Pursuing dynamics of minimal residual leukemic subclones in relapsed and refractory acute myeloid leukemia during conventional therapy.

BACKGROUND: Acute myeloid leukemia (AML) is characterized by clonal heterogeneity, leading to frequent relapses and drug resistance despite intensive clinical therapy. Although AML's clonal architecture has been addressed in many studies, practical monitoring of dynamic changes in those subclones during relapse and treatment is still understudied. METHOD: Fifteen longitudinal bone marrow (BM) samples were collected from three relapsed and refractory (R/R) AML patients. Using droplet digital polymerase chain reaction (ddPCR), the frequencies of patient's leukemic variants were assessed in seven cell populations that were isolated from each BM sample based on cellular phenotypes. By quantifying mutant clones at the diagnosis, remission, and relapse stages, the distribution of AML subclones was sequentially monitored. RESULTS: Minimal residual (MR) leukemic subclones exhibit heterogeneous distribution among BM cell populations, including mature leukocyte populations. During AML progression, these subclones undergo active phenotypic transitions and repopulate into distinct cell population regardless of normal hematopoiesis hierarchic order. Of these, MR subclones in progenitor populations of patient BM predominantly carry MR leukemic properties, leading to more robust expansion and stubborn persistence than those in mature populations. Moreover, a minor subset of MR leukemic subclones could be sustained at an extremely low frequency without clonal expansion during relapse. CONCLUSIONS: In this study, we observed treatment persistent MR leukemic subclones and their phenotypic changes during the treatment process of R/R AML patients. This underscores the importance of preemptive inhibition of clonal promiscuity in R/R AML, proposing a practical method for monitoring AML MR subclones.

Prediction of metabolites associated with somatic mutations in cancers by using genome-scale metabolic models and mutation data.

BACKGROUND: Oncometabolites, often generated as a result of a gene mutation, show pro-oncogenic function when abnormally accumulated in cancer cells. Identification of such mutation-associated metabolites will facilitate developing treatment strategies for cancers, but is challenging due to the large number of metabolites in a cell and the presence of multiple genes associated with cancer development. RESULTS: Here we report the development of a computational workflow that predicts metabolite-gene-pathway sets. Metabolite-gene-pathway sets present metabolites and metabolic pathways significantly associated with specific somatic mutations in cancers. The computational workflow uses both cancer patient-specific genome-scale metabolic models (GEMs) and mutation data to generate metabolite-gene-pathway sets. A GEM is a computational model that predicts reaction fluxes at a genome scale and can be constructed in a cell-specific manner by using omics data. The computational workflow is first validated by comparing the resulting metabolite-gene pairs with multi-omics data (i.e., mutation data, RNA-seq data, and metabolome data) from acute myeloid leukemia and renal cell carcinoma samples collected in this study. The computational workflow is further validated by evaluating the metabolite-gene-pathway sets predicted for 18 cancer types, by using RNA-seq data publicly available, in comparison with the reported studies. Therapeutic potential of the resulting metabolite-gene-pathway sets is also discussed. CONCLUSIONS: Validation of the metabolite-gene-pathway set-predicting computational workflow indicates that a decent number of metabolites and metabolic pathways appear to be significantly associated with specific somatic mutations. The computational workflow and the resulting metabolite-gene-pathway sets will help identify novel oncometabolites and also suggest cancer treatment strategies.

Analysis of clinical and genomic profiles of therapy-related myeloid neoplasm in Korea.

BACKGROUND: Therapy-related myeloid neoplasm (T-MN) rarely occurs among cancer survivors, and was characterized by poor prognosis. T-MN has germline predisposition in a considerable proportion. Here, clinical characteristics and germline/somatic variant profiles in T-MN patients were investigated, and the findings were compared with those of previous studies. METHODS: A review of medical records, cytogenetic study, targeted sequencing by next-generation sequencing, and survival analysis were performed on 53 patients with T-MN at a single institution in Korea. RESULTS: The patients were relatively younger compared to T-MN patients in other studies. Our T-MN patients showed a high frequency of complex karyotypes, -5/del(5q), and -7/del(7q), which was similar to the Japanese study group but higher than the Australian study group. The most common primary disease was non-Hodgkin lymphoma, followed by breast cancer. The detailed distributions of primary diseases were different across study groups. Seven patients (13.2%) harbored deleterious presumed/potential germline variants in cancer predisposition genes (CPG) such as BRIP1, CEBPA, DDX41, FANCM, NBN, NF1, and RUNX1. In the somatic variant profile, TP53 was the most frequently mutated gene, which was consistent with the previous studies about T-MN. However, the somatic variant frequency in our study group was lower than in other studies. Adverse factors for overall survival were male sex, older age, history of previous radiotherapy, previous longer cytotoxic therapy, and -5/del(5q). CONCLUSION: The findings of our study corroborate important information about T-MN patients. As well as a considerable predisposition to CPG, the clinical characteristics and somatic variant profile showed distinctive patterns. Germline variant testing should be recommended for T-MN patients. If the T-MN patients harbor pathogenic germline variants, the family members for stem cell donation should be screened for carrier status through germline variant testing to avoid donor-derived myeloid neoplasm. For the prediction of the prognosis in T-MN patients, sex, age, past treatment history, and cytogenetic findings can be considered.

Interpretation of EBV infection in pan-cancer genome considering viral life cycle: LiEB (Life cycle of Epstein-Barr virus).

We report a novel transcriptomic analysis workflow called LiEB (Life cycle of Epstein-Barr virus) to characterize distributions of oncogenic virus, Epstein-Barr virus (EBV) infection in human tumors. We analyzed 851 The Cancer Genome Atlas whole-transcriptome sequencing (WTS) data to investigate EBV infection by life cycle information using three-step LiEB workflow: 1) characterize virus infection generally; 2) align transcriptome sequences against a hybrid human-EBV genome, and 3) quantify EBV gene expression. Our results agreed with EBV infection status of public cell line data. Analysis in stomach adenocarcinoma identified EBV-positive cases involving PIK3CA mutations and/or CDKN2A silencing with biologically more determination, compared to previous reports. In this study, we found that a small number of colorectal adenocarcinoma cases involved with EBV lytic gene expression. Expression of EBV lytic genes was also observed in 3% of external colon cancer cohort upon WTS analysis. Gene set enrichment analysis showed elevated expression of genes related to E2F targeting and interferon-gamma responses in EBV-associated tumors. Finally, we suggest that interpretation of EBV life cycle is essential when analyzing its infection in tumors, and LiEB provides high capability of detecting EBV-positive tumors. Observation of EBV lytic gene expression in a subset of colon cancers warrants further research.

A Case of Pulmonary Langerhans Cell Sarcoma Simultaneously Diagnosed with Cutaneous Langerhans Cell Histiocytosis Studied by Whole-Exome Sequencing.

Langerhans cell histiocytosis (LCH) and Langerhans cell sarcoma (LCS) are clonal proliferations of Langerhans-type cells. Unlike in LCH, the pathophysiology and clinical course of LCS are unclear due to its rarity. Here, we report the case of a 73-year-old male patient who was diagnosed with cutaneous LCH and pulmonary LCS at the same time. Pathological review of these 2 tumors revealed similar immunohistochemical findings. However, the tumor cells in LCS had more aggressive cytological features than those in LCH. Results of BRAF mutation analysis using real-time PCR were negative for both tumors. In whole-exome sequencing (WES), stop-gain mutations in TP53 gene were discovered only in LCS cells. The mechanism of development of LCS from various progenitor cells is currently unclear. According to the results of the WES study, changes in TP53 gene might have contributed to the malignant features of LCS.

Biophysical characterization of soluble Pseudomonas syringae ice nucleation protein InaZ fragments.

Ice nucleation protein (INP) with its functional domain consisting of multiple 48-residue repeat units effectively induces super-cooled water into ice. Circular dichroism and infrared deconvolution analyses on a soluble 240-residue fragment of Pseudomonas syringae InaZ (InaZ240) containing five 48-residue repeat units indicated that it is mostly composed of ß-sheet and random coil. Analytical ultracentrifugation suggested that InaZ240 behaves as a monomer of an elongated ellipsoid. However, InaZ240 showed only minimum ice binding compared to anti-freeze proteins. Other P. syringae InaZ proteins with more 48-residue repeat units were made, in which the largest soluble fragment obtainable was an InaZ with twelve 48-residue repeat units. Size-exclusion chromatography analyses further suggested that the overall shape of the expressed InaZ fragments is pH-dependent, which becomes compact as the numbers of 48-residue repeat unit increase.

Tracking solvent and protein movement during CO2 release in carbonic anhydrase II crystals.

Carbonic anhydrases are mostly zinc metalloenzymes that catalyze the reversible hydration/dehydration of CO2/HCO3 (-) Previously, the X-ray crystal structures of CO2-bound holo (zinc-bound) and apo (zinc-free) human carbonic anhydrase IIs (hCA IIs) were captured at high resolution. Here, we present sequential timeframe structures of holo- [T = 0 s (CO2-bound), 50 s, 3 min, 10 min, 25 min, and 1 h] and apo-hCA IIs [T = 0 s, 50 s, 3 min, and 10 min] during the "slow" release of CO2 Two active site waters, WDW (deep water) and WDW' (this study), replace the vacated space created on CO2 release, and another water, WI (intermediate water), is seen to translocate to the proton wire position W1. In addition, on the rim of the active site pocket, a water W2' (this study), in close proximity to residue His64 and W2, gradually exits the active site, whereas His64 concurrently rotates from pointing away ("out") to pointing toward ("in") active site rotameric conformation. This study provides for the first time, to our knowledge, structural "snapshots" of hCA II intermediate states during the formation of the His64-mediated proton wire that is induced as CO2 is released. Comparison of the holo- and apo-hCA II structures shows that the solvent network rearrangements require the presence of the zinc ion.

Crystallization and preliminary X-ray crystallographic analysis of carboxyl-terminal region 4 of SigR from Streptomyces coelicolor A3(2).

Full-length SigR from Streptomyces coelicolor A3(2) was overexpressed in Escherichia coli, purified and submitted to crystallization trials using either polyethylene glycol 3350 or 4000 as a precipitant. X-ray diffraction data were collected to 2.60 Šresolution under cryoconditions using synchrotron X-rays. The crystal packs in space group P43212, with unit-cell parameters a=b=42.14, c=102.02 Å. According to the Matthews coefficient, the crystal asymmetric unit cannot contain the full-length protein. Molecular replacement with the known structures of region 2 and region 4 as independent search models indicates that the crystal contains only the -35 element-binding carboxyl-terminal region 4 of full-length SigR. Mass-spectrometric analysis of the harvested crystal confirms this, suggesting a crystal volume per protein weight (VM) of 2.24 Å3 Da(-1) and 45.1% solvent content.

The crystal structure of an activated Thermotoga maritima CheY with N-terminal region of FliM.

In bacterial chemotaxis, the levels of phosphorylated CheY in association with FliM determine the sense of the flagella rotation, which in turn controls the bacterial swimming behavior. We report the 1.7Å resolution crystallographic structure of the Thermotoga maritima BeF(3)(-)-activated CheY in complex with the CheY-binding N-terminal region of FliM. Analysis of the structure in comparison to the previously reported Escherichia coli counterpart reveals that similar regions of H4-ß5-H5 in CheY and the helix in FliM are used for the complex interfaces. Our structure also indicates that the correlated movement of Phe101 and Ser82 (F-S coupling) in T. maritima CheY upon phosphorylation and FliM binding, parallels that of Tyr106 and Thr87 (Y-T coupling) demonstrated in E. coli CheY. Furthermore, significant displacements of the ß4-H4 loop in both CheYs impose a crucial role of this loop, which can be related to flagellar switch component binding or to propagating changes that is necessary during the CheY-mediated reversal of the motor.